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LONG-TERM  GOALS 

The  long-term  goals  of  the  study  are  to  model  and  measure  bedforms  in  tidal  inlets  and  river  mouths. 
Specifically,  the  goals  are:  1)  to  use  an  existing  self-organization  model  to  predict  multiple  scales  of 
bedform  formation,  growth  and  migration  in  combined  steady,  tidal  and  wave-driven  flows,  2)  to 
compare  model  predictions  with  measurements  (from  both  the  literature  and  from  a  tidal  inlet/river 
mouth  experiment),  3)  to  test  the  hypothesis  that  subaqueous  bedforms  will  grow  indefinitely  for  a 
given  set  of  conditions  and  that  different  scales  of  bedforms  will  occur  simultaneously  because  of  this 
continuous  growth,  and  4)  to  incorporate  this  bedform  model  into  community  modeling  systems  (eg 
CSDMS)  to  improve  modeling  of  flow,  sediment  transport  and  morphology  change  in  river  mouths 
and  tidal  inlets. 

OBJECTIVES 

The  specific  objectives  of  this  study  have  been  to 

>  develop  and  adapt  the  present  model  for  flows  in  river  mouths  and  tidal  inlets  (expand  the 
model  to  2-D  flows,  improve  flow  drivers,  scale  the  model  up  for  larger  spatial  domains  and 
begin  to  examine  multiple  bedform  scales). 

>  continue  working  with  the  CSDMS  so  that  the  present  model  can  be  utilized  by  that  community 
modeling  environment. 

>  participate  in  the  ONR-supported  tidal  inlet  and  river  mouth  experiment  with  Tom  Lippmann  at 
CCOM/UNFI  to  generate  a  data  set  on  bedforms  in  a  tidal  inlet. 

APPROACH 

Bedforms  in  sandy  environments  are  ubiquitous,  occurring  in  rivers,  river  mouths,  estuaries,  tidal  inlets 
and  on  open-coast  beaches.  Bedforms  act  as  roughness  elements,  altering  the  flow  and  creating 
feedback  between  the  bed  and  the  flow.  In  doing  so,  they  are  intimately  tied  to  erosion,  transport  and 
deposition  of  sediments  (eg  Parsons  et  al.  2005,  Emstsen  et  al.  2005).  It  has  been  suggested  that 
bedforms  in  rivers  and  tidal  inlets  are  dynamically  similar  to  Aeolian  dunes  and  bedforms  on  the 
continental  shelf  and  in  the  surf  zone  (Best  2005,  Frank  and  Kocurek  1996,  Nemeth  et  al.  2007, 
Gallagher  2003).  Because  of  this  similarity,  Gallagher  (2011)  developed  a  model  for  bedforms  in  the 
nearshore,  based  on  the  principles  of  work  by  Werner  (1995),  who  hypothesized  that  Aeolian  dunes 
were  self-organized  features  and  as  such  could  be  modeled  with  a  relatively  simple  model. 


It  has  been  suggested  that  self-organization  is  responsible  for  the  formation  of  many  different  types  of 
morphological  patterns,  including  river  meanders  (Stolum  1996),  sorted-patterned  ground  (Kessler  & 
Werner  2003),  beach  cusps  (Coco  et  al.  2000),  wind  ripples  (Nishimori  &  Ouchi  1993)  and  Aeolian 
dunes  (Wemer  1995,  Reffet  et  al.  2010).  In  each  of  these  pattern-forming  systems,  complexity  arises 
from  nonlinear  interactions  between  the  system  and  the  environment,  from  dissipative  processes  such 
as  friction,  turbulence  and  sediment  transport,  and  from  being  open  (both  material  and  energy  are 
exchanged  across  system  boundaries)  and  therefore  never  in  equilibrium  (Wemer  1999). 

Werner  (1995)  used  a  ‘hierarchical’  approach  (Ahl  &  Allen  1996)  to  modeling  self-organized  systems, 
wherein  processes  at  different  temporal  and  spatial  scales  are  distinct  from  each  other  and  can  be 
separated.  With  this  approach,  grain-scale  sediment  transport  is  parameterized  with  simple  rules  to 
drive  bedfonn-scale  dynamics.  Gallagher  (2011)  developed  a  hierarchical  model  to  predict  nearshore, 
combined  flow  megaripples.  The  model  consists  of  a  matrix  of  sediment  slabs  that  represent  a  spatial 
domain  or  a  region  of  a  bed  across  which  sediment  is  moving.  The  sand  slabs  are  picked  up  and  moved 
according  to  a  transport  model  (either  simple  rules  similar  to  Werner  (1995)  or  a  physics-based 
formulation,  e.g.  Bailard  1981,  Ribberink  1998).  Sediment  transport  is  driven  by  the  free  stream 
velocity,  u,  which  is  modeled  with  a  sinusoidal  velocity,  a  measured  velocity  signal  from  the  natural 
surf  zone  or  with  a  Rayleigh  distributed  wave  velocity  signal.  At  each  time  step,  the  flow  is  the  same  at 
all  locations  in  the  domain  except  for  an  imposed  random  spatial  fluctuation  representing  local 
turbulence.  However,  once  bedforms  are  created,  the  flow  around  the  bedforms  is  altered  via  feedback: 
flow  is  reduced  in  the  lee  of  a  bedfonn  to  simulate  a  velocity  shadow  zone  and  flow  is  accelerated  over 
the  crest  of  a  large  bedfonn.  These  spatial  alterations  to  the  flow  generate  gradients  in  transport,  which 
alter  the  bed.  Feedback  is  required  for  bedfonn  growth  and  development  (Gallagher  2011).  In 
addition,  the  slope  of  the  bed  is  not  allowed  to  exceed  17°. 

The  long-term  plan  for  this  research  is  to  work  with  and  adapt  the  self-organization  model,  developed 
for  nearshore  bedforms.  This  model  is  being  adapted  for  predicting  bedforms  in  the  combined  flows  of 
tidal  inlets  and  river  mouths.  In  these  environments,  oscillatory  flows  with  wave  frequencies  are 
superimposed  on  the  quasi-steady  flows  associated  with  tides  (oscillatory  but  with  a  much  longer 
period  than  the  surf  waves)  as  well  as  steady  flows  (possibly  with  seasonal  variations)  exiting  river 
mouths.  These  complex,  but  naturally  realistic,  flows  are  being  incorporated  to  predict  the  growth  and 
migration  of  dunes  and  the  evolution  of  multiple  scales  of  bedforms.  In  addition  to  combined  flows 
and  multiple  scales,  variations  owing  to  spatially  varying  grain  size  will  also  be  examined.  This  model 
lends  itself  to  tackling  these  dynamically  complex  issues,  because  relatively  simple  changes  can  be 
implemented  to  test  the  importance  of  factors  such  as  lateral  flows,  feedback  changes,  grain  size  and 
subtle  3-D  morphology  changes.  Model  results  will  be  compared  with  data  from  the  literature  (eg, 
Hanes  2012  ,  Jerolmack  and  Mohrig  2005,  Ernstsen  et  al.  2005)  and  with  data  collected  as  part  of  the 
River  Mouths  and  Tidal  Inlets  DRI  experiments  in  collaboration  with  Tom  Lippmann  (UNH),  Steve 
Elgar  (WHOI)  and  Peter  Traykovski  (WHOI). 

Lastly,  the  model  has  been  submitted  to  the  Community  Surface  Dynamics  Modeling  System.  CSDMS 
is  a  community  of  experts  promoting  open-source  modeling  of  earth  surface  processes.  They  develop, 
support,  and  disseminate  integrated  software  modules  that  predict  the  movement  of  fluids  and  the  flux 
of  sediments  and  solutes  in  landscapes  and  sedimentary  basins  (for  more  information  see  their  website 
at  http://csdms.colorado.edu/wikEIntroduction).  Here,  the  bedfonn  model  can  be  used  and  improved 
by  others  and  integrated  into  larger-scale  fluid  and  morphodynamic  models  with  the  intention  of 
improving  predictions  of,  bed  roughness,  wave  and  flow  dissipation,  and  sediment  transport. 


WORK  COMPLETED 


Significant  progress  has  been  made  in  adapting  the  bedfonn  model  for  river  mouths  and  inlets  (see  Fig 
1  for  a  schematic).  The  flow  model,  which  drives  sediment  transport  and  bedfonn  dynamics,  has  been 
extended  from  1-D  to  2-D.  Preliminary  results  suggest  that  this  advance  has  improved  model 
predictions  significantly  by  adding  complexity  and  slowing  the  bedfonn  formation  process  (Fig  4). 
Other  improvements  to  the  model  flow  field  that  have  also  improved  predictions  are  the  addition  of  a 
Rayleigh  distributed  wave  field  and  the  addition  of  flow-dependent  turbulence. 

The  original  MATLAB  version  of  the  model  is  available  on  the  CSDMS  website.  The  model  has  been 
translated  from  MATLAB  to  Fortran  to  fascilitate  increasing  the  domain  size  and  running  on  the 
CSDMS  High  Performance  Computing  Cluster.  The  model  is  now  being  run  successfully  (and 
speedily)  for  larger  domains.  In  addition,  with  a  Fortran  version,  the  model  can  be  included  as  a 
‘component’  in  CSDMS  making  it  easily  integratable  with  other  models.  Currently,  I  am  sitting  in  on  a 
computer  science  course  to  learn  the  Python  programming  language.  The  CSDMS  staff  scientists  have 
recommended  Python  as  a  more  modem  and  flexible  language  for  modeling.  Thus  I  am  moving  in  the 
direction  of  translating  the  model  into  Python  as  well.  (I  am  learning  a  lot!)  These  steps  are  considered 
neccesary  in  the  model  development,  although  in  themselves  they  do  not  produce  interesting  or 
measureable  results. 


Figure  1.  Flowchart  of  model  improvement  completed  and  underway. 

In  addition  to  these  coding/computer  science  aspects  of  the  model  improvement  a  number  of  model 
adaptations  are  being  considered  to  improve  the  model’s  accuracy.  With  the  increased  domain  size 
now  possible,  the  model  needs  to  have  spatially  varying  flow  and  bed  domains  to  realistically 
reproduce  the  larger  scale  tidal  inlet  regions  to  be  modeled.  (In  contrast  to  the  smaller  nearshore 
megaripple  domain,  where  the  flow  was  spatially  uniform.)  In  addition,  the  model  feedback  is  being 
adapted  to  create  more  realistically  shaped  bedfonns  and  to  (hopefully)  improve  the  predicted  growth 


dynamics  (Fig  4).  In  this  case,  a  crude  mechanism  for  suspended  load  bypassing  will  be  included  in  the 
feedback/transport  routines  (this  is  in  its  infancy). 

I  presented  the  model  at  the  ICCE  2012,  AGU  Fall  Meeting  2012  and  received  lots  of  excellent 
feedback.  I  hope  to  establish  a  collaboration  with  modelers  at  UPC  (Universitat  Politecnica  de 
Catalunya.  We  have  discussed  a  meeting  next  summer  during  the  RCEM  conference).  The  model  was 
also  presented  at  the  ONR  program  review  in  September  2012  in  Denver  and  at  the  New  River  Inlet 
Experiment  (see  below)  Review  in  April  2013  (despite  no  funding). 

I  participated  in  the  tidal  inlet  experiment  at  New  River  Inlet  with  Tom  Lippmann  in  May  2012. 
Preliminary  bedform  data  from  that  experiment  (collected  using  a  multibeam  sonar  to  make  detailed 
bed  measurements  and  using  rotary  side  scan  (Traykovski))  are  exciting  and  suggest  multiple  scales  of 
bedforms  exist,  that  bedforms  vary  significantly  depending  on  depth,  sediment  availability  and  that 
these  bedforms  are  changing  constantly  as  the  waves  and  tides  change.  Some  examples  of 
Traykovski’s  preliminary  data  can  be  seen  at  http://vimeo.com/44806773.  In  these  movies,  multiple 
scales  of  bedforms  are  observed  to  form,  change  directions,  be  destroyed  and  reform  over  the  course  of 
a  tidal  cycle.  Through  collaboration,  these  data  sets  will  be  used  to  test  the  model. 


Predicted  bedforms  using  the  self-organization  model 
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Figure  2.  Top  row:  modeled  bedforms  for  combined  (oscillatory  and  steady)  flows  (left),  oscillatory  flow  only  (middle)  and 
steady  flow  only  (right).  Oscillatory  flows  have  a  period  of  10  secs  and  an  amplitude  of  75  cm/s  (left)  and  95  cm/s  (middle). 
The  steady  flows  were  20  cm/s  (left)  and  50  cm/s  (right).  Bottom  row:  observed  bedforms.  Red  arrows  show  the 
correspondence  between  similar  modeled  and  observed  features. 


RESULTS 

This  model  was  developed  for  nearshore  flows  (combined  waves  and  currents)  and  sediment  transport. 
Predicted  nearshore  bedforms  are  shown  in  Fig  2  (top  left)  and  are  similar  to  observed  features  in  the 
surf  zone  (Fig  2,  bottom  left,  from  Clarke  and  Werner  2004.  See  also  Gallagher  2003).  Altering  the 
forcing  flow  field  in  a  model  of  this  type  to  simulate  other  flows  (eg,  quasi-steady  tidal  flow,  steady 
fluvial  flow  or  purely  oscillatory  flow)  is  simple.  For  example,  by  setting  the  steady  flow  to  Om/s  and 
running  the  model  with  purely  oscillatory  flow,  the  combined-flow  forms  lose  any  directionality  (eg, 
bedform  ends  no  longer  point  downstream)  and  become  irregularly  shaped  lumps  (Fig  2  top  middle). 
This  change  in  shape  is  also  illustrated  in  the  profiles  of  bedforms  (Fig  3)  taken  along  the  solid  lines  in 


Fig  2.  The  combined  flow  features  show  some  directionality,  with  steep  downstream  faces  and 
broader,  more  shallowly  sloped  upstream  slopes  (Fig  3  a),  but  when  the  steady  flow  is  removed  they  are 
no  longer  asymmetric  (Fig  3b).  Fluvial  bedforms  are  simulated  by  setting  the  oscillatory  amplitude  to  0 
m/s  and  using  only  steady  flow.  The  predicted  steady  flow  bedforms  (Fig  2  top  right)  are  highly 
asymmetric  (Fig  3c)  and  similar  in  planform  to  observed  bedforms  (Fig  2,  bottom  right). 


Distance  (m) 

Figure  3.  Profiles  across  simulated  bedforms  along  lines  marked  in  Fig  2(top  panels).  In  c)  steady  flow  pro  files  are  highly 
asymmetric,  while  in  a)  that  asymmetry’  is  reduced  by  a  superimposed  oscillatory  current.  In  b),  with  no  steady  flow  at  all, 

bedforms  are  symmetric  in  profile. 


The  model  predicts  that  bedforms  begin  as  random  irregularities  on  the  bed  and,  via  feedback  between 
the  flow  and  the  bed,  coalesce  into  small  bedforms.  As  bedforms  continue  to  evolve,  smaller,  faster 
bedforms  merge  with  larger,  slower  ones,  causing  crest-  and  wave-lengths  to  grow.  Thus,  younger 
bedforms  tend  to  be  short-crested,  shorter  in  wavelength  and  irregular  in  shape,  while  more  mature 
features  are  longer  in  both  wavelength  and  crest  length.  This  merging  and  lengthening  is  observed  in 
nature  (eg,  Clarke  and  Werner  2004)  and  in  other  modeling  studies  (eg,  Coco  and  Murray  2007, 
Werner  and  Kocurek  1999,  Jerolmack  and  Mohrig  2005).  Clarke  and  Wemer  (2004)  observed  that  the 
growth  of  bedforms  (wavelengths)  from  a  flat  bed  to  maturity  occurred  first  linearly  and  then  became 
logarithmic.  This  corresponds  to  the  theoretical  model  of  Wemer  and  Kocurek  (1999),  which  attributed 
the  dynamics  and  growth  of  bedforms  to  the  behavior  of  defects  or  the  ends  of  bedform  crests.  The 
change  in  growth  (from  linear  to  logarithmic)  was  attributed  to  the  growth  and  lengthening  of  bedform 
crests  and  the  reduction  in  numbers  of  defects.  The  present  model  predicts  this  transition  in  growth  rate 
(and  the  reduction  in  defect  number),  suggesting  that  the  model  captures  well  the  bedform  dynamics. 
The  natural  megaripples  in  Clarke  and  Werner  (2004)  made  the  transition  at  around  12  hrs.  The 
modeled  megaripples  grow  and  transition  much  more  quickly:  25  mins  for  sinusoidal  flows  and  50 
mins  for  natural  measured  velocities  (Fig  4). 

The  difference  in  growth  rate  of  modeled  megaripples  driven  with  the  sinusoidal  versus  the  natural 
velocities  may  be  explained  by  examining  the  velocity  records.  The  largest  amplitudes  of  the  natural 
cross-shore  velocity  from  the  measured  time  series  are  over  100  cm/s  and  the  root  mean  square  (RMS) 
is  32  cm/s.  The  sinusoidal  flows  have  amplitudes  of  75  cm/s  and  a  RMS  of  65  cm/s.  This  difference  is 
because  the  measured  velocities  are  skewed  (with  the  strongest  flows  having  a  short  duration)  and 
irregular,  with  the  largest  velocities  (>75  cm/s)  occurring  infrequently.  So,  under  natural  flows,  high 
transport  rates  are  intermittent.  In  contrast,  the  sinusoidal  flows  reach  their  maximum  velocity  every 
cycle  and  drive  high  rates  of  sand  transport  consistently.  Therefore,  bedforms  are  built  more  quickly 
under  the  consistent  sinusoidal  flows  and  more  slowly  under  the  variable  natural  flows.  Neither  flow 


field  reproduces  the  natural  growth  rate  and  transition  time  of  12  hrs  observed  by  Clarke  and  Werner 
(2004).  The  long  transition  time  observed  in  the  natural  surf  zone  likely  results  from  the  even  higher 
variability  of  the  total  flow  field,  including  more  realistic  turbulence,  more  realistic  acceleration  on  the 
bedform  crest  (acting  to  reduce  amplitude  growth),  variation  in  direction,  variation  in  tidal  level,  which 
Clarke  and  Werner  (2004)  state  is  the  dominant  controller  of  the  magnitude  of  the  depth-dependent, 
wave-driven  flows  in  the  surf  zone  (Raubenheimer  2002),  and  possibly  the  frequent  interruption  of  the 
feedback  mechanisms  by  turbulence  from  breaking  waves. 


Figure  4.  Bedforms  wavelength  growth  with  time.  The  black  arrow  marks  the  time  for  which  Clarke  and  Werner  (2004) 
found  nearshore  megaripples  transitioned  from  early  linear  growth  to  later  logarithmic  growth  (12  hrs).  The  red  and  blue 
arrows  mark  the  transition  times  for  sinusoidally  (25  mins)  and  natural-flow  (50  mins)  driven  bedforms,  respectively.  The 
magenta  points  are  for  bedforms  driven  by  a  directionally  varying  flow  with  the  new  2D  model.  There  is  no  quantitative 
transition  time  for  these  preliminary  data,  but  it  can  be  seen  that  they  are  growing  more  slowly  than  their  ID  counterparts. 


As  part  of  the  present  study,  adaptations  have  been  made  to  the  model  to  include  more  realistic  flows 
and  to  allow  the  model  to  be  scaled  up  (Fig  1).  For  example,  by  implementing  2-D  flows  and  using  the 
y-direction  measured  flows  that  correspond  to  the  natural  x-direction  flows  used  in  the  calculations  for 
Fig  4,  the  model  was  run  and  the  model  growth  rate  was  slowed  further  (the  magenta  points  fall  below 
and  to  the  right  of  the  blue  points  in  Fig  4).  Note  that  this  model  run  has  not  been  completed  for  the 
long  time  needed  to  quantify  this  result  so  the  transition  time  has  not  yet  been  determined.  These 
preliminary  results  support  the  idea  that  variability  of  the  bedform-building  flows  are  important  in 
understanding  the  growth  and  dynamics  of  the  bedforms. 

IMPACT/APPLICATION 

This  model  is  being  adapted  and  applied  to  different  environments.  At  this  time  the  model  is  being 
expanded  to  be  able  to  predict  bedforms  in  larger  scale  more  highly  variable  environments.  For 
example,  the  flow  field  and  conditions  varyed  significantly  in  time  and  in  space  during  the  New  River 
Inlet  experiment,  including  the  deep  sediment-starved  channel  with  very  strong,  quasy  steady,  tidal 
flows  and  the  waves-plus-tidal-flow-dominated  shallow,  sandy  shoal  that  included  breaking  and 
directionally  varying  flows.  Soon,  the  model  will  be  compared  with  observations  from  these  varied 
environments.  This  will  be  the  first  attempt  at  modeling  tidal  inlet  and  river  mouth  bedforms  with  the 
self  organization  model.  It  is  expected  that  a  simple  model  of  this  type  could  be  expanded  to  model 


other  environments.  By  beginning  to  work  with  CSDMS  at  this  time,  it  is  hoped  that  this  model  will  be 
easily  integrated  into  larger-scale  flow  and  morphology  models  and  will  help  improve  the  predictive 
capabilities  of  hydro-  and  morpho-dynamics  in  general. 
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